Wave Propagation in Overhead Wires 
with Ground Return 



By JOHN R. CARSON 



I V HE problem of wave propagation along a transmission system 
■*■ composed of an overhead wire parallel to the (plane) surface of 
the earth, in spite of its great technical importance, does not appear 
to have been satisfactorily solved. 1 While a complete solution of 
the actual problem is impossible, on account of the inequalities in the 
earth's surface and its lack of conductive homogeneity, the solution 
of the problem, where the actual earth is replaced by a plane homo- 
geneous semi-infinite solid, is of considerable theoretical and practical 
interest. The solution of this problem is given in the present paper, 
together with formulas for calculating inductive disturbances in 
neighboring transmission systems. 

The axis of the wire is taken parallel to the s-axis at height h 
above the xz-p\ane and passes through the y-axis at point 0' as 
shown in Fig. 1 herewith. The "image" of the wire is designated 
by 0". 

For y>0 (in the dielectric) the medium is supposed to have zero 
conductivity, while for -y<0 (in the ground) the conductivity of the 
medium is designated by X. The xs-plane represents the surface of 
separation between dielectric and ground. 

We consider a wave propagated along the s-axis and the current, 
charge and field are supposed to contain the common factor 
exp ( — Tz-\-iwt), which, however, will be omitted for convenience in 
the formulas. The propagation constant, T, is to be determined. It is 
assumed, ab initio, as a very small quantity in c.g.s. units. 2 

In the ground (y^O) the axial electric force is formulated as the 

1 See Rudenberg, Zt. f. Angewandt, Math. u. Mechanik, Band 5, 1925. In that 
paper the current density in the ground is assumed to be distributed with radial 
symmetry. The resulting formulas are not in agreement with those of the present 
paper. Since this paper was set up in type I have learned that formulas equivalent 
to equations (26), (28), (31) for the mutual impedance of two parellel wires were 
obtained by my colleague, Dr. \\. A. Campbell, in 1917. It is to be hoped that 
his solution will be published shortly. 

1 The simplifying assumptions introduced in this analysis are essentially the same 
as those employed and discussed in "Wave Propagation Over Parallel Wires: The 
Proximity Effect," Phil. Mag., Vol. xli, April, 1921. 
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general solution, symmetrical with respect to x, of the wave equation; 
thus 



E.. 



= - f F{u) cos xfi fl+Wdii, y^O 



(1) 



where 



a = 47rXco, 

X= conductivity of ground in elm. c.g.s. units, 



i=V-l, 
w/27r = frequency in cycles per second. 

(In the following analysis and formulas, elm. c.g.s. units are employed 
throughout). 
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Assuming that in the ground E x and E y are negligible compared 
with E 3 , we have from the formula, curl E= —-H, 

iull x = — — E z 
dy 



i(j}Hv = —-E z . 

dx 
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Whence, in the ground 

H x =i-f vV+«« .F(n). cos xii^F+todn, ( 2 ) 

//v=i f fi.F(n). sm x/x.e^^+^dfi, (3) 

it being understood that y^O. The function F(n) in the preceding 
formulas is to be determined by the boundary conditions. 

In the dielectric, II X and II y may be regarded as composed of two 
components; thus 

H X =H X »+H X ', 

where II X °, II y ° designate the field due to the current I in the wire, and 
IIx , Hy' the field of the ground current. 

Neglecting axial displacement currents in the dielectric, and as- 
suming that the wire is of sufficiently small radius so that the distri- 
bution of current over its cross section is symmetrical, we have 

H<0= cos_e' 2J> 

p 

ff, - 5 ^. 21, (5) 

p 

where p = \/x 2 + (y — h)~, 

cos9' = ^=2 (0) 

P 

sin Q' = x/p. 
The secondary magnetic field H x , II y ' is taken as 

Xoo 
<j>((i)cos xp.e-wdp t (7) 

0(/x)sin xn.e-wdn. (8) 

At the surface of separation y=0, II X °, H y ° are expressible as the 
Fourier integrals 

H X ° = 2I f cos Xfi.e-todp, (9) 

IIy = 21 f sin x/i.e-^dn. (10) 
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Also at the surface of separation of the two media (y = 0), H x and II y 
must be continuous. Equating the values of H x and II y at y = 0, 
as given by (2), (3) and by (7), (8) and (9), (10), we have 

}-H.F{tx)=2I.e-^-<i>^), 
whence 

W = ;===— 4/, (11) 

0(m) = 5^+^2)r*".2/, (12) 

VV + ia + m 

which determines the functions F(n) and #(/■»)• 

Inserting the value of F(n), as given by (11) in (1), the axial electric 
force E- in the ground (v^O) and therefore the distribution of current 
density in the ground is expressed as a Fourier integral in terms of the 
frequency a>/2ir, the current / in the wire, the height h of the wire 
above ground, and the conductivity X of the ground. Similarly the 
insertion of 4>(m)> as given by (12) in formulas (7) and (8) gives the 
magnetic field II X , II y in the dielectric. Thus 

/°° e-vk 

/ = ; . , Wjn+to cos Xfx.dn, y^O. (13) 

This can be further simplified if we write 

x' =X\/a 

y' =y\/a 

h' = h\/'a, 
whence 

£ s =-4a>/^ {y/f+i-tie-'^./^ 1 ^ cosx'txdy., y^O. (14) 

The axial electric force in the dielectric is now to be formulated. 
This is always derivable from a vector and a scalar potential; thus 

E z =-io>A z -^-V, (15) 

where A z is the vector potential of the axial currents, and V the 
scalar potential. Consequently, 



(18) 
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E z (x,y)-E z (x,0) = - i»(M*&) -M*,0)) -|^( V(x,y) -V ). (16) 

Here E-(x,0) is the axial electric intensity at the surface of the ground 
plane (y = 0), and 

A s (x,y) -A s (x,0) = P H x (x,y)dy. (17) 

V(x,y) — Vo is the difference in the scalar potential between the point 
x,y and the ground, which is due to the charges on the wire and on 
the surface of the ground. For convenience, it will be designated by V. 
By means of (16) and the preceding formulas we get 3 

(Vr+i-n)e~ lh ' +y ' )}l cos *V dn 

- i 2u>I log (p"/ P ') ~V, y^O 

where 

P ' = V(h- y y-+x n - 

= distance of point x,y from wire, 

p'^Vih+yy+x* 

= distance of point x,y from image of wire. 

The first two terms on the right hand side of (18) represent the 

electric force due to the varying magnetic field; the term V 

dz 

represents the axial electric intensity due to the charges on the surface 
of the wire and the ground. If Q be the charge per unit length, V 
is calculable by usual electrostatic methods on the assumption that 
the surface of the wire and the surface of the ground are equipotential 
surfaces, and their difference of potential is QIC where C is the electro- 
static capacity between wire and ground. 4 

II 

By aid of the preceding analysis and formulas, we are now in a 
position to derive the propagation constant, T, and characteristic 
impedance, K, which characterize wave propagation along the system. 
Let z denote the "internal" or "intrinsic" impedance of the wire per 

3 As a check on this formula note that together with (14) it satisfies the condition 
of continuity of E, at y =0. 

4 See "Wave Propagation Over Parallel Wires: The Proximity- Effect," Phil. Mag., 
Vol. xli, Apr., 1921. 
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unit length. (With small error this may usually be taken as the resist- 
ance per unit length of the wire.) The axial electric intensity at the 
surface of the wire is then zl. Equating this to the axial electric in- 
tensity at the surface of the wire as given by (18) and replacing 3/3z 
by — r, we have 

zl = -4W f™ Wf+i-tie-^dn 

J ° (19) 

-t 2wJ.log (p" /a) +TV. 

Writing V=Q/C and 

iaQ = rl-GV=TI-^Q, 

where G is the leakage conductance to ground per unit length, we 
have, solving for T, 

r 2 = (G+ioC)[z+i2u. log ( P "/a)]+4co f" (V?+7-|i)«^ w *«)- ( 2 0) 

Writing this in the usual form 

r- = (R+iX)(G+ia>C), (21) 

the characteristic impedance is given by 

X 2 = ^±^ (22) 

and the series impedance per unit length of the circuit is 

R+iX = Z = z+i2o. log (p"/a)+4<o C {\/J+l-n)e-* h, Hn. (23) 

It will be observed that the first two terms on the right hand side of 
(23) represent the series impedance of the circuit if the ground is a per- 
fect conductor; the infinite integral formulates the effect of the finite 
conductivity of the ground. 

The mutual impedance* Z i2 between two parallel ground return 
circuits with wires at heights hi and hi above ground and a separation 
x between their vertical planes is given by 

z u -a». log (p"/p')+4«X~ (VT+i-tie-* 1 ™'' cos *'* *■» W 

6 It will be noted that the mutual impedance is equal to the axial electric intensity 
at the axis of the second wire due to the varying magnetic field of unit current in 
the first wire and its accompanying distribution of ground current. 
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where 

p" = V(/*i+/* 2 ) 2 +* 2 
p'=V(hi-h 2 y+x* 

h% =Jhy/ a. 
x' = x\/~a. 



From the preceding the series self impedance of the ground return 
circuit may be conveniently written as 

Z = Z°+Z' (25) 

and the mutual impedance as 

Zi2 = Zi2-\-Z'i2 (26) 

where Z°, Z° ]2 are the self and mutual impedances respectively, on the 
assumption of a perfectly conducting ground, and 

Z' = ±o>f™ {y/W+i-tie-^dv, (27) 

Zj 2 = 4co f°° (Vl^+i-u)e~ lhl ' +h - ) ' i cos x'a da. (28) 

•A) 

The calculation of the circuit constants and the electromagnetic 
field in the dielectric depends, therefore, on the evaluation of an 
infinite integral of the form 

J(p,q)=J= / (Vu 2 -\-i-u)e-"f l . cos qu du. (29) 

In terms of this integral 

Z' = 4w.7(2/i'.0) (30) 

Z , i2 = 4co./(/*/+/* 2 , 1 x'). (31) 

To the solution of the infinite integral J(p,q) we now proceed. 

Ill 

The solution of equation (29), that is, the evaluation of J{p,q) 
can be made to depend on the solution of the infinite integral 



/ y/ifi+c?.e-todi 



546 BELL SYSTEM TECHNICAL JOURNAL 

which has been worked out and communicated to me by R. M. Foster. 
It is 



^{^ 1 (a^) + G(a/3)[ 



where K\{x) is the Bessel function of the second kind and first order 
as defined by Jahnke und Emde, Funktionentafeln, pg. 93, and 
G(x) is the absolutely convergent series 

# 2 x* x 6 



3 2 .5 ' 3 2 .5 2 .7 

On the basis of this solution, it is a straightforward though intricate 
and tedious process to derive the following solution for J(p,q) of 
equation (29). 

Writing r=Vp 2 +q 2 and e = tan~ 1 (q/p), it is J = P+iQ 

in which 

P = | (1 - 54 ) + l(,ogA), 2+ |e.,/ 



1 _L 1 4. X 



(32) 



^ = T+T( l0 ^f,) (1 - 54) -T e ^ 
i 1 » _i 1 x 

V2 8 V2 * 



(33) 



In these equations log 7 is Euler's constant: 

7 = 1.7811, log —=0.11593, log 7 = 0.57722 and a,, <r«, <x 3 , <r h s 3l ft', 
7 

5 4 , s/, are infinite series defined as follows: 

*-i^t)~»-s1bGj) ,00,, » + -- ■■ 
^=ik!(T) %in2e -3ir!(T) r ' sin6e+ ---' 
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1 (r \ 4 1 (r \ 8 

• V4 ' = 2!3A2 J sin4G "4!5!U J sin 8e + • • • • 



r cos r 5 cos 5G , r° cos 9G 

o\ = 



(T3 = 



3 3 2 .5 2 .7 ' 3 2 .5 2 .7 2 .9 2 .11 ' 

r 3 cos 39 r 1 cos 70 , r 11 cos 119 



3 2 .5 3 2 .5 2 .7 2 .9 3 2 .5 2 .7 2 .9 2 .11 2 .13 " ' ' 
ff -( I +5-4 I )li ! (0 2cos2e 

-( 1 +l+3- + 5-?)3Hi(5)' cosfle+ --' 



5 

= jSo approximately, 

. -( i+ Ws + l +1 s-^)Md' cosSe+ ■■■■ 

5 
— ^s^ approximately. 

It is to be regretted that the foregoing formulas appear so com- 
plicated. The series, however, are very rapidly convergent and for 
r^2 only the two leading terms of each series need be retained. For 
r^l, only the leading terms are of importance. 

For the important range r^l/4, 

P = l jrr cos 9 + ^ cos 29(o.6728+log-)+^9sin 29, (34) 

o 3v 2 1(J \ r ' lo 

Q= -0.0386+1 log (y) + ~r cos 9. (35) 

For r>5 the following asymptotic expansions, derivable from (29) 
by repeated partial integrations, are to be employed. 

_1 cos 9 _ cos 29 1 cos 39 3 cos 59 .„ . 

~" -y/2 r r 2 y/2 r 3 y/% r b ' ' ' ' 

r>5 

n _l cos 9 1 cos 39 3 cos 59 .„ . 
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For large values of r(r>10), these reduce to 
l+»cos6 cos 26 



V2 



-,r>10. 



(38) 



In view of the somewhat complicated character of the function in 
the range l/4^r^5 the curves shown below have been computed. 

These show J = P+iQ as a function of r for = 0, — , — , — , — . By 

interpolation it is possible to estimate with fair accuracy the value 
of the functions for intermediate values of 6. 
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Fig. 2 P= real part of J 
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Fig. 3 Q = imaginary part of J 



IV 



The theory and formulas of the preceding sections will now be 
reviewed and summarized as regards their principal applications to 
technical transmission problems where the ground forms, in whole 
or part, the "return" part of the circuit. In such problems we are 
interested in the electric intensity in the dielectric and in the ground, 
and in particular in the self impedance and mutual impedances of 
ground return circuits. The calculation of these quantities is pro- 
vided for by the general analysis and the solution of the infinite 
integral J. Reference should be made to Fig. 1 shown in section I 
for the geometry of the system and coordinate system employed. 

1. The Axial Electric Intensity E z in the Dielectric. (See equations 
(15) and (18) ). 



E z = -% V-(i2a, log 0>"/p')+W)7 
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where 

p' = V(h-y)*+x* 
= distance of point in dielectric from axis of wire. 

P " = V(h+y) 2 +x* 

= distance of point in dielectric from, image of wire. 

r = P "Va 

e = sin- i (x/ P ") 

a=47rXco. 

These values of r and Q are, of course, to be employed in calculating 
J = P-\-iQ from the formulas and curves of the preceding section. 
As a special case the electric intensity at the surface of the earth is 

E 3 =-4coJI 

p"=\//?b ! 

r = P "V« 

e = sin- l (x/p") 

2. Self Impedance of Ground Return Circuit. (See equations (25) . 

(27), (30) ). 

Z = Z°+4co/ 
Z" = self impedance with perfectly conducting ground. 
r = 2h \/a 

e = o. 

3. Mutual Impedance of Ground Return Circuits. (See equations (26), 

(28), (31) ). 

Z 12 = Z^+4co/ 

Z°2 = mutual impedance with perfectly conducting ground. 



r- xfaV(h 1 +h 2 ) 2 +x?=p"V<x 

6 = sin- 1 (x/p"). 

The axial electric intensity E z in the ground (y <0) is given by equation 
(1), and the subsequent analysis, whence 

E z = -4co/ f "(V?+?-|i).cob x'n.e-W+'^'^dii 
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where, as before 

x' = x y/a 

h' = hy/'a 
i nd 

g'= -y/a times the depth below the surface of the ground. 

= gVoc. 

The integral can undoubtedly be evaluated in somewhat the same way 
as (29) and can in any case be numerically computed without much 
difficulty. Owing, however, to the secondary technical interest in the 
electric intensity below the surface of the earth, the detailed solution 
has not been undertaken, nor has the magnetic field been worked out. 

V 

The practical utility of the preceding theory and formulas will now 
be illustrated by a brief sketch of their application to two important 
transmission problems. 

The Wave Antenna 

When a transmission line with "ground return" is employed as a 
radio receiving antenna it is called a wave antenna. The theory and 
design of such an antenna requires a knowledge of the transmission 
characteristics of the ground return circuit, which are calculable, as 
shown above, from the geometry and constants of the overhead wire, 
together with Z' = icoJ, which may be termed the "ground return" 
impedance. 

We assume that the wire is approximately 30 ft. above the ground 
(h=10 3 ) and that the frequency is 5.10 4 c.p.s. corresponding to the 
frequency employed in Trans-Atlantic radio communication. The 
ground conductivity X is exceedingly variable, depending on the 
locality and weather conditions. Calculations of Z' will therefore 
be made for two extreme cases, X = 10~ 12 and X = 10 -14 which should 
cover the range of variation encountered in practice. 

For X = 10- 12 , 
y/a= \/47rX^ = 2. 10" 3 
and for X = 10- 14 , 
a = 2.l0- 4 . 
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Correspondingly, r = 2h\/a has the values 4.0 and 0.4, respectively. 
Reference to the preceding formulas and curves for J, for r = 4.0 and 
r = 0A, give 

7=0.126+*' 0.168, X = 10~ 12 

J = 0.323+*' 0.871, X = 10~ 14 

whence the corresponding values of Z r are 

Z' = 4co. (0.126+* 0.168), 

Z' = 4co. (0.323+* 0.871). 

These are the "ground return" impedances per unit length in elm. 
c.g.s. units; to convert to ohms per mile they are to be multiplied by the 
factor 1.61 X10 -4 . Consequently setting co = 7r.l0 5 , we get 

Z' = 6.44tt(1.3+*'1.7), X = 10- 12 

Z' = 6.44tt (3.2+* 8.7), X = 10-' 4 . 

Comparison of these formulas shows that an hundred-fold increase 
in the resistivity of the ground increases the resistance component of 
the ground return impedance by the factor 2.5 and increases its 
reactance only five-fold. That is to say, the ground return impedance 
is not sensitive to wide variations in the resistivity of the earth, a 
fortunate circumstance in view of its wide variability and our lack of 
precise information regarding it. 

Induction from Electric Railway Systems 

A particularly important application of the preceding analysis is 
to the problems connected with the disturbances induced in parallel 
communication lines by alternating current electric railways. As- 
suming the frequency as 25 c.p.s., we have corresponding to X = 10 -12 
and X = 10- 14 , 

^= 0.45 X 10" 4 and 0.45 X 10" 5 . 



Taking the height of the trolley wire as approximately 30 ft., h = 10 3 
and assuming the parallel telephone as the same height above ground 
and separated by approximately 120 ft., # = 4.10 3 , and 

r=v / aV(2A) 2 -^ 2 
= 4.47X10 3 V^ 
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and corresponding to the values of a taken above 

r = 0.2 and 0.02 in round numbers, 
while 

4 
G = sin l — y= = 63° 30' approximately. 

For both cases, therefore, we can employ, in calculating J = P-\-iQ, 
the approximate formulas, 



P = ^"^ r - cos0+ ^ cos2e (- 6728+log 7) + 5 esin 



29, 



Q= -0.0386+ \ log (-) H X - F -r cos 9. 

For X = 10 - ' 2 and r = 0.2, this gives 

7 = 0.369+*' 1.135 
and 

Zi 2 =4w(0.369+* 1.135). 

The foregoing assumes that the only return conductor is the ground. 
If, however, an equal and opposite current flows in the rail we must 
subtract from the foregoing mutual impedance, the mutual impedance 
between rail and telephone line; that is, the mutual impedance Z' 2 
between the telephone line and a conductor at the surface of the 
earth. For this case 

P "= V/^+^ 2 = 4.12Xl0 3 
9 = sin- 1 ^^ = 76° 

cos 9 = 0.242, r = 0.184 for X=10-' 2 . 
The corresponding value of J is 
7 = 0.378+*" 1.165 

and the resultant mutual impedance between railway and parallel 
telephone line is, 

Z , 12 = 4o J (0.369-0.378+i (1.135-1.165) ) 
= 4co(0.009-z 0.030). 

The very large reduction in mutual impedance, due to the current 
in the rail, is striking. 
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For the case of X = 10 -14 , the corresponding calculations give 
Z{ 2 =4w(0.391+*2.27) 
with no current in rail, and 

Zi2=4w( - 0.001 -* 0.002) 

with equal and opposite current in rail. It is evident from these 
figures that the reduction in mutual impedance, due to the current 
in the rail, is practically independent of the ground conductivity, at 
least at the separation specified. 



